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In two dimensions a microscopic theory providing a basis for the naive analogy between a quantized 
vortex in a superfluid and an electron in a uniform magnetic field is presented. Following the 
variational approach developed by Peierls, Yoccoz, and Thouless, the cyclotron motion of a vortex 
is described by the many-body wave function, which is a linear combination of Feynman wave 
functions centered at different positions. An integral equation for the weighting functions of the 
superposition is derived by minimizing the energy functional. The matrix elements of the kernel are 
the overlaps between any two displaced Feynman wave functions. A numerical study is conducted 
for a bosonic superfluid based on a Hartree ground state. A one-to-one correspondence between 
the rotational states of a vortex in a cylinder and the cyclotron states of an electron in the central 
gauge is found. Like the Landau levels of an electron, the energy levels of a vortex are highly 
degenerate. However, the gap between two adjacent energy levels does not only depend on the 
quantized circulation, but also increases with energy, and scales with the size of the vortex. The 
fluid density is finite at the vortex axis and the vorticity is distributed in the core region. The 
effective mass of a quantized vortex defined by the inverse of the energy-level spacing is shown to 
be logarithmically divergent with the size of the vortex. 

PACS numbers: 67.40.Vs, 02.30. Rz 



I. INTRODUCTION 

Quantized vortices have been known to play a signif- 
icant part in the behaviors of a supetfluid for several 
decades since OnsagerQ and Feynman.0 Understanding 
the vortex dynamics from the quantum theoretical point 
of view still remains as a challenging problem after these 
years. Recently, an analogy with the cyclotron motion 
of an electron has helped us to make a few progresses on 
this subject. P 

The phenomenological theoryjj which models a quan- 
tized vortex in a two-dimensional superfluid as a charged 
particle in a gauge field, has been very intriguing for un- 
derstanding the dynamics of the vortex. This analogy is 
motivated by the fact that, in classical physics, the mo- 
tion of a vortex in an ideal fluid is similar to the motion 
of an electron in an uniform magnetic field. In both cases 
the broken time reversal symmetry leads to a transverse 
force proportional to the velocity. The Magnus force due 
to the Bernoulli pressure difference is equal to the mass 
density of the fluid times the vector product between the 
circulation vector along the vortex axis and the velocity 
of the vortex relative to the background fluid. The well- 
known Lorentz force is equal to the charge of the electron 
times the vector product between the magnetic field and 
the velocity of the electron relative to the field. Follow- 
ing the simple relationship between the energy and the 
momentum of an ideal—fluid, a vortex is like a particle 
with an effective mass.Ea The contributions to the effec- 
tive mass come both from the particles trapped inside 
the vortex core and from the fluid disturbed by the mo- 
tion of the vortex. As a result of the velocity-dependent 
transverse force, the trajectory of a classical vortex con- 
sists of a circular orbit around a guiding center with an 



arbitrary radius but a fixed angular frequency. In the su- 
perfluid 4 He, the dynamics of a quantized vortex is also 
governed by the Magnus force and the effective mass be- 
cause the compressibility and the viscosity of the fluid 
are comparably small 

Having the strong analogy in classical physics between 
an electron and a vortex in mind, one expects to find 
some common features in their fully quantized theories. 
However, the actual theoretical developments of these 
two systems are quite different from one another. On 
one hand, enormous progress has been made on the quan- 
tum theory of electrons in the presence of a high mag=. 
netic field since the discovery of the quantum Hall effect. El 
Even though the real system consists of strongly corre- 
lated electrons, the basic physics of the integer quantum 
Hall effect can be understood from the simple theory of a 
non-interacting electron gas on a two-dimensional plane. 
In order to compare with the vortex motion later, I will 
discuss the solution in the central gauge so that the angu- 
lar momentum in the perpendicular direction to the plane 
is conserved. In the central gauge, electrons move as if 
in a rotating frame with an effective simple-harmonic po- 
tential. The energy eigenstates shift up and down from 
the simple harmonic spectrum according to their indi- 
vidual angular momentum, and form highly degenerate 
Landau levels. The detailed solution to this idealized 
problem without any disorder is given in the appendix 
of this paper. On the other hand, there is no quantum 
theory starting from first principles which supports the 
existence of the cyclotron motion and degenerate levels 
of vortex states. Most theories of vortex dynamics are 
semi-classical in the sense that the motion of a vortex 
is described explicitly by the coordinates of the vortex 
center moving on a classical trajectory. The governing 
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equations for the vortex center and the quantization con- 
ditions are not known due to a fundamental difference 
between an electron and a vortex: an electron is a real 
particle obeying Schrodinger equation, but a vortex is a 
collective object in a many-body system. Experimental 
evidences for these collective modes in the superfiuid 4 He 
have been shown in a more complicated situation, where 
the cyclotron motion is coupled to a traveling wave along 
the vortex axis.Q In high-T c superconductors, the vortex 
cyclotron resonance, has also been studied in several re- 
cent experiments. Em 

It is the purpose of the present paper to provide a 
microscopic theory for this phenomenological analogy, 
and show that the rotational states of a quantized vor- 
tex also form highly degenerate energy levels similar to 
the Landau levels in the integer quantum Hall effect. 
In order to describe these rotational states, which are 
collective modes in a many-body system, I construct 
variational wave functions using the projection method 
adapted from thp_.nuclear theory. It was suggested by 
Hill and WhceleiO that the wave functions of rotational 
states in a nucleus could be constructed by integrating 
out the surface variables, which specify the position and 
the orientation of the nucleus. The wave functions of a 
quasi-harmonic oscillator were proposed for the weighting 
funetioas of the integration. Peierls, Yoccoz and Thou- 
lessOlla later pointed out that the weighting functions 
could be determined from the variational principle. I 
apply their scheme on a vortex in a bosonic superfiuid, 
and calculate the energy spectrum and the weighting 
functions in the dilute-gas limit. In certain sense these 
weighting functions can be understood as the "single- 
particle wave functions" of the vortex. At the end of 
this paper, the weighting functions are shown numeri- 
cally to be identical to the eigenfunctions of a rotating 
two-dimensional harmonic oscillator, which is consistent 
with the analogy of cyclotron motion. Degenerate energy 
levels like Landau levels are also found in this theory. 
However, the energy gap between two neighboring levels 
is not a constant, and scales logarithmically with the size 
of the vortex. 

In Sec. H a brief review is given on the semi-classical 
wave functions of a quantized vortex in a superfiuid. 
These many-body wave functions are constructed with 
the help of classical hydrodynamics, and will be used as 
the basis functions for constructing the new trial wave 
functions. 
In Sec. 
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the projection method is discussed in terms 
of the quantum fluctuation of the vortex position. I pro- 
pose the new set of trial wave functions based on the 
linear combination of the semi-classical wave functions. 
An integral equation for the weighting function of the su- 
perposition is obtained by using the variational principle. 

In Sec. IV a system of bosons with short-ranged re- 
pulsive interactions is considered. In the dilute-gas limit, 
the ground state of the system can be described by the 
Hartree approximation. The kernel of the integral equa- 
tion is reduced to a sum of products of one-particle inte- 



grals. 

In Sec. |v| the exact solutions of the dilute system is ob- 
tained by solving the integral equation numerically. En- 
ergy levels formed by lots of degenerate states are found. 
A partial explanation for the degeneracy is given. To vi- 
sualize the structure of a vortex, I calculate the radial 
distributions of the number density and the current den- 
sity. Since a simplified ground state is used, the curves 
only show qualitative features comparing; to the results of 
quantum Monte Carlo calculations ,Ej~li3 which, however, 
are only avai lable for a static vortex. 

In Sec. VI an excellent agreement between the vortex 
motion and the cyclotron motion is found by comparing 
the weighting functions and the electron wave functions. 
In analogy to the cyclotron motion, the effective mass of 
a vortex is link to the inverse of the energy gap, which 
scales logarithmically with the size of the vortex. 



II. SEMI-CLASSICAL WAVE FUNCTIONS 

At the zero temperature, the wave function for an uni- 
form flow in a superfiuid can be obtained by performing 
a Galilean transformation on the ground state, 
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where m is the mass of the particle, v s is the flow ve- 
locity, Yj is the position vector of each particle, N is the 
total number of particles, and is the ground state 
wave function. The ground state wave function is cho- 
sen to he real and positive apart from an arbitrary phase 
factor £3 The wave function in Eq. (0) can be generalized 
to represent a flow slowly varying in space, 
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where the velocity field is described by the gradient of 
the phase of the wave function, 



v s (r) = -V0(r) . 
m 



(3) 



Consequently, the motion of the superfiuid has to be a 
potential flow in a simply connected region. 

In a rotating superfiuid, vortices with quantized circu- 
lations are formed. A naive picture of a vortex is a core 
with jjon- vanishing vorticity immersed in an irrotational 
nuidia For a single static vortex fixed at a point ro in an 
infinite homogeneous system, the velocity field outside 
the core is inversely proportional to the distance, 
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where 0(r — ro) is the unit vector in the azimuthal direc- 
tion relative to the point ro. In this region a wave func- 
tion of the form in Eq. (0) can be constructed with the 
phase function <f>{r) equal to the azimuthal angle 6(r— ro). 
Inside the core, this construction scheme fails completely 
because the velocity field is no longer curl-free, and can- 
not be incorporated into a phase function with a single 
variable. Moreover, the actual structure of the core is in 
general not known. In the superfluid 4 He, the size of the 
core is about the size of a helium atom or less. Needless to 
say, it is very difficult to probe the core structure at this 
scale experimentally. Theoretically, the simplest remedy 
for constructing a trial wave function is to assume that 
the phase function possess the same form everywhere, but 
the amplitude of the wave function is reduced around the 
core, 
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where g(r) is a cut-off function, which goes to zero at 
the origin, and goes to one aswnptotically. Eq. (|5|) 
was first discussed by Feynman, a and was extensively 
used in literatures with the cut-off function determined 
variationallyJla This wave function is semi-classical in the 
sense that the vortex is exactly at the point ro, like a 
(5-function in quantum mechanics. Nevertheless, it is a 
fairly good description for a vortex in the superfluid 4 He 
because the uncertainty in the vortex position is on the 
order of the inter-atomic spacing. 

In the situation where the vortex is moving with a ve- 
locity Vo, there is an additional dipolar flow generated by 
the motion of the core. The exact flow pattern depends 
on the detail of the core. For example, the velocity po- 
tential for a vortex with a hard cylindrical core is 
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where a is the radius of the hard core.Q The additional 
dipolar flow given by the second term in the above equa- 
tion is called backflow. The backflow is quite important 
for considering the dynamical properties of the vortex. 
In classical hydrodynamics the backflow is responsible 
for the effective mass of a vortex. One can compute the 
total kinetic energy of the fluid from the velocity field. 
The first-order term in vo vanishes, and the coefficient 
of the second-order term gives the effective mass, which 
in two dimensions is equal to the amount of the fluid 
displaced by the hard core, 



M c s = pwa , 
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where p is the mass density of the fluid. In the super- 
fluid 4 He, the backflow needs to be included in the trial 
wave functions in order to obtain a quantitative agree- 
ment between calculations and experimental data on the 
energy spectrum of quasi-particle excitations, especially 



on the part of rotonsS The rotons are usually visualized 
as vortex rings. 

In the rest of this paper I have completely neglected 
all the backflow corrections in the trial wave functions for 
the technical difficulty. As a result, only qualitative fea- 
tures of the energy spectrurtt-afld of the effective mass can 
be obtained. A few paperstJ'Ej have included the back- 
flow in the trial wave function, but only the static vortex 
is discussed there. In the static case, the correction is 
small because the vortex is not really moving besides the 
zero-point motion. 



III. PROJECTION METHOD 

In the Feynman wave function a vortex is described 
both by the coordinates of particles and by the coordi- 
nates of the vortex core. An intuitive way to consider the 
vortex motion is to make the coordinates of the core time- 
dependent like real dynamical variables. However, like 
the coordinates for the center of mass, the coordinates 
of the core are collective variables which do not show up 
explicitly in the underlying microscopic Hamiltonian. A 
set of dynamical variables including both the particle co- 
ordinates and the core coordinates is a redundant set. In 
the absence of external pinning centers, the motion of the 
vortex must have been accounted in a proper quantum 
mechanical description of the iV-particle system. How 
does the motion of the vortex appear as part of the mo- 
tion of a superfluid rather than as something imposed 
from outside? Conversely, how can the motion of the vor- 
tex be described by degrees of freedom which are not real 
dynamical variables? Similar questions were asked long 
ago in a different context regarding the rotational states 
of a large deformed nucleus. A large nucleus is a many- 
body system which under certain conditions can move or 
rotate as a whole. In the framework of the Hartree-Fock 
theory, the wave function of a nucleus is written in terms 
of a Slater determinant of one-particle wave functions in 
a self-consistent potential well. If the translation or the 
rotation of the nucleus is considered through the position 
or the orientation of the potential well, then one also has 
a redundant set of dynamical variables. 

The answer to the above questions was first given by 
Hill and Wheeler in terms of the quantum fluctuation 
of the potential well. In the context of vortex dynam- 
ics, the corresponding answer is the fluctuation in the 
vortex position. In a collective motion, the particle vari- 
ables and the collective variables are strongly correlated. 
Since there are quantum fluctuations in the motion of the 
particles in a superfluid, there must be fluctuations in the 
position of the vortex, too. Mathematically, we integrate 
out the fluctuating variables other than the particle co- 
ordinates, and obtain a linear combination of Feynman 
wave functions with given core coordinates, 
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where the weighting function G(ro) represents the quan- 
tum fluctuation. This wave function depends only on the 
coordinates of particles, and in the mean time provides 
a description of the vortex motion through the weight- 
ing function. {i"o} is called generator coordinates, and 
is a good set of independent coordinates in the Hilbcrt 
space because any two Feynman wave functions centered 
at different positions are nearly orthogonal to each other. 
In this case, a naive guess for the weighting function is 
the wave function of a rotating two-dimensional harmonic 
oscillator. I will show later that this guess is indeed the 
correct one. 

An alternative point of view to Eq. (^) was later pro- 
vided by Peierls, Yoccoz and Thouless, who also devel- 
oped a method to determine the weighting function sys- 
tematically. A vortex located at a certain point in space 
breaks the translational symmetry explicitly. There- 
fore we have a family of semi-classical wave functions 
parametrized by their center coordinates, and all of them 
give the same expectation value in energy. Generally 
speaking, the variational energy can be further reduced 
by using a superposition of these degenerate functions as 
the new trial wave function, and the weighting function 
can be determined by minimizing the energy integral, 



G(r ) = f n ,i{r )e 



no,, 



(12) 



(E) = 



J/d 2 r^ 2 r G*(r-)G(r )(^|g|vI/ r 



(9) 



where H is the Hamiltonian. The trial wave function 
can also regain the proper symmetry of the Hamiltonian 
with a suitable choice of the weighting function. This 
approach is then called the projection method because it 
projects out a state with the desired symmetry. 

I start with a system of N bosons in a disk geometry 
with radius R. The generic Hamiltonian is 

h 2 N N 

H = E t v * + + 2 £ r ^ < ( 10 ) 

i—l — \ 

where the particles interact with one another through 
the potential V(r), and the boundary of the system is 
represented by the confining potential Vb(r), which is a 
step function with an infinite height. The explicit form of 
the inter-particle potential V(r) is not directly relevant to 
our discussion here, which will become clear in a moment. 

The first step is to find the velocity potential of a single 
vortex located at ro = (xo,yo)- This problem is solved 
by using an image vortex with the opposite-circulation 
located at a distance R 2 /ro from the origin,E!l 
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One can easily verify that Eq. ( pT| ) is the solution by 
checking that the radial current vanishes on the bound- 
ary. Because of the cylindrical symmetry of the system, 
the weighting function in Eq. (||) must be of the form, 



where n is the principal quantum number, and / is the 
angular momentum. Following the convention of Feyn- 
man,Eil I can rewrite the trial wave function in the fol- 
lowing form, 
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where all the phase factors are combined together, 
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For simplicity, I have set the cut-off function g(r) in 
Eq. (||) to be unity. In this case, even though the Feyn- 
man wave function ^S ro becomes singular, the trial func- 
tion in Eq. (j^) remains to be regular. The cut-off function 
cannot be set to unity if one wants to include the back- 
flow in the trial wave function. The expectation value of 
the Hamiltonian relative to the ground state energy E gs 
is 
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It is a generic feature for trial wave functions of the form 
in Eq. (13|) that the inter-particle potential does not show 
up explicitly in the variational energy e n> i. The influ- 
ences from the inter-particle interactions only come in 
indirectly through the ground-state wave function. By 
varying the weighting function f n ,i{f'o) to minimize £ n ,U 
I obtain an integral equation for fn t i(ro), 
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where the kernel of this integral equation contains two 
parts: J; is the angular average of the overlap between 
any two Feynman wave functions located at the radius 
ro and r' , and K\ is the similar overlap weighted by the 
kinetic energy, 
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where the primed functions $J and come from the 
complex conjugate of the trial wave function, and are 
referred to the corresponding unprimed functions with 
ro substituted by r . It is clear that the integrand only 
depends on the relative angle, A6q = 8' — 6q. All particle 
coordinates in the integrand are equivalent as a result of 
the Bose symmetry. 
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where O is the opening angle spanned from r — r' to 
r — ro , and 0/ is the corresponding angle with respect to 
the image vortices. In this form, one can easily see that 
the integrand is independent of the coordinate system. 
These angles can be easily formulated by considering the 
coordinates as complex variables, 



IV. HARTREE APPROXIMATION 

In order to simplify the kernel further and obtain an 
analytic form, one needs an explicit form of the ground- 
state wave function. In the Hartree approximation for 
bosons, the ground-state wave function is broken down 
to a product of normalized one-particle wave functions, 



N 



(20) 



where u(r) satisfies the non-linear Schrodinger equation 
with a self-consistent potential, 

V sc (r) = (N - 1) J d 2 r'V(\r - r'|Mr') 2 + V h (r) . (21) 

For a dilute system of bosons with short-ranged repulsive 
interactions, this one-particle wave function is approxi- 
mately constant except for a transition layer near the 
boundary. Assuming that the effective repulsive inter- 
actions is infinitely strong, I can neglect the transition 
layer, and set my ground state to a constant. Since the 
higher-order correlations among particles are completely 
neglected in the Hartree approximation, the A^-particle 
multiple integral in Eqs. ( jig ) and (19) is reduced to a 
simple product of one-particle integrals, 
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where the normalized ground state is A~ N / 2 , and A is 
the area of the system. 

From now on, I will set the unit length to be the inter- 
particle spacing, a = R/y/N, which is the only length 
scale in this model. Other length scales can only come in 
implicitly through the spatial variations of the ground- 
state wave function, and are discarded in my simple ap- 
proximation. For a two-dimensional helium film at the 
saturation density, a is about 2.7A from numerical calcu- 
lations.^ The corresponding energy scale associated with 
this length scale is h 2 /ma 2 ss 1.7 K. 

The phase factor in the one-particle integrals in Eqs. 
( p2| ) and (^3|) can be rewritten in terms of the opening 
angles between the two generator coordinates, r and r , 
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One begins the calculation with the assumption that the 
system size R is sufficiently large so that r /R can be 
treated as a small parameter. This assumption is self- 
consistent only if the weighting functions are localized 
far away from the boundary. The validity of this as- 
sumption can be checked later as we discuss the solu- 
tions. The common overlap integral in both Ji and Ki 
is evaluated by dividing the integration range into three 
distinct regions. The first region is the circle centered at 
r c = (ro + r )/2 with diameter d — \tq — r |, and has 
an area A\ = nd 2 /A. In this region, one can regard the 
phase contributed by the image vortices as a constant, 
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The remaining part of the integral in the large R limit is 
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Since the integrand is dimensionless, the integral must 
scale as d 2 , and the numerical coefficient is calculated 
explicitly in Appendix B. The second region is the ring 
centered at r c from the inner radius d/2 to the outer 
radius R — r c , and has an area Ai = ir(R — r c ) 2 — A\. 
In this region, it is convenient to shift the origin of my 
coordinate system to r c . The variables in Eqs. (Eq) and 
(^6| ) will be changed accordingly. I can expand e with 
respect to (ro/r) and (r' Q /r), and expand e lQ/ with re- 
spect to (rro/R 2 ) and (rr' /R 2 ) up to the second order. 
Keeping only the terms, which survive after the angular 
integration, I obtain 
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For the radial integral, the upper bound is taken to be 
R, the lower bound is d/2, and the integral gives 
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where the factor 5/2 is contributed by the image vortices. 
The last region is simply the remaining area A3 = A — 
A\ — A2 like a thin crescent moon. In this region, one 
can again regard the phase contributed by the images as 
a constant, and the remaining integral gives 
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In the second line of the above equation, I have changed 
the integration range from A3 to A so that the cross terms 
linear in r vanish. After putting Eqs. (p8|), ( |3C| ) and j3l| ) 
together, the original integral is obtained, 
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where a = 7/8 + 1.6639/4tt = 1.0074. The overall phase 
factor contributed by the image vortices is closely related 
to the fact that lh is the total angular momentum of 
the system, not the angular momentum of the effective 
motion of the vortex. This connection will become clear 
as we discuss the solutions in the next section. 

The other overlap integral in Eq. (|2^) is weighted by 
the kinetic energy. The kinetic energy term can also be 
rewritten in terms of the opening angles, 
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In this case the image vortices only contribute to the 
overall phase factor. I can again shift the origin to r c , 
and the first approximation to the integral is 



r 2 * r R dr 
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where e l9 w cosG w 1. The correction to this result 
must be a dimensionless number since there is no other 
dimensional parameters other than d, 
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where (3 — —0.30685. The detailed calculation of this 
number can be found in Appendix B by evaluating the 
integral numerically. Therefore, we reach the final form 
of the kernel in the large R limit as, 
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This kernel is sharply peaked when the two generator 
coordinates are close to one another because both the 
separation d and the phase of Eq. ( ^8|) are small. Notice 
a peculiar feature that the kernel depends on the system 
size explicitly. 



V. SOLUTIONS OF THE INTEGRAL EQUATION 

Even though I have greatly simplified the system, the 
kernel still remains in an integral form. Neither the an- 
alytic form of the kernel, nor the series solutions of the 
integral equation have been worked out because the last 
integral contains logarithms. Nevertheless, the integral 
equation can be solved numerically, which leads to some 
interesting results. In fact, solving this homogeneous in- 
tegral equation is basically like solving the eigen-system 
of a matrix. The variational energy e n j, and the weight- 
ing function f n j are the eigenvalues and eigenvectors of 
the matrix J^ l Ki. The generator coordinates, ro and r' , 
now only take discrete values, and act like the matrix in- 
dices. The only complication is that each matrix clement 
of J\ or Ki is a one-variable integral. 

Instead of working with the full-size matrix, in which 
ro is running from to R, it is sufficient to work only 
with a sub-matrix, in which ro is running only from 
to a smaller cut-off R c , for two reasons: (a) The ma- 
trix elements far away from the diagonal are consider- 
ably smaller than the elements around the diagonal. This 
statement is equivalent to say that two far apart Feynman 
wave functions are essentially orthogonal to each other, 
(b) Since the kernel is sharply peaked, it is reasonable 
to assume that the weighting functions are localized in 
the region where the kernel is peaked. The use of such 
cut-off is self-consistent for the eigenvectors, which are 
localized in the region bounded by the cut-off. This self- 
consistent condition is the same as the one for calculating 
the kernel Eqs. (^) and (|37|). The greatest advantage of 
this cut-off is that one can maintain the same numeri- 
cal accuracy for calculations on fairly large system sizes 
without drastically increasing the sizes of the matrices. 
A typical matrix in my calculation has about two hun- 
dred elements on each side for R c — 10. The calculation 
is carried out for a series of different sizes of matrices at 
a fixed i? c , and an extrapolation to the continuum limit 
is performed in the end to obtain the final results. 

Numerical results are summarized in Fig. [j] and Fig. |^. 
Fig. [l] represents a few examples of the weighting func- 
tions, whose modulus squares are normalized to unity. 
These weighting functions are well-localized, so that the 
self-consistent condition is justified. Their functional 
forms are found to be independent of the system size 
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as long as R is large enough to satisfy the self-consistent 
condition. In general, the weighting functions should de- 
pend both on the system size and the angular momen- 
tum, because the kernel depends on both parameters in a 
non-trivial way. Besides the evidence from numerical so- 
lutions, I am unable to provide any analytical calculation 
here to prove this independence rigorously. Physically, 
this independence means that the motion of a vortex is 
not sensitive to the boundary conditions, which is quit 
reasonable. Fig. |^ is the energy spectrum for various 
different system sizes. The angular momentum index is 
dropped because the energy eigenvalues in the continuum 
limit are independent of the angular momentum. These 
large degeneracies in the variational energies can be un- 
derstood approximately, in the large R limit, as follows. 
For clarity, I will only show the equations for J; because 
the equations for K\ are almost identical. For weighting 
functions sufficiently far away from the origin, the ker- 
nel of their integral equation is peaked in a small region 
around r « r' w y/N — I with \N - l\ > 1. The kernel 
can be approximated by 



Ji 



dye 



-(In x /2R/d+a)d 2 p -iy[r' -(N-l)/r' ] 



(39) 



where the variables are changed to x = \r' Q — ro|, y = 
ra Ado an d d 2 w x 2 + y 2 . Two kernels with different an- 
gular momenta are related to each other in the following 

way, 



Ji+Ai(r' + Ar ; ,r ) w Ji(r' ,r ) 



(40) 



For a given Al, one can find that An = y/N — I — Al — 
y/N — I by keeping the phase of the integrand in Eq. j39| ) 
invariant when r is near y/N — I. Therefore, there are 
approximate relations among the variational energies and 
among the weighting functions, 



£n,l+Al ~ 7 

fn,l+Ai(ro + Ar ) w f n ,i(r ) 



(41) 
(42) 



As a result, the variational energies are highly degener- 
ate, and the degeneracies are of order N. The similarity 
between any two weighting functions at the same energy 
level with different angular momenta is demonstrated by 
comparing the plots with I = N — 9 and with I = N — 10 
in Fig. 0. 

To see the structure of a vortex in this model, one 
can compute the number density and the current density. 
The number density is 



N 



i=i 



*n,i« (43) 



Since the density is independent of the azimuthal an- 
gle by symmetry, I can perform an angular average, and 
rewrite the integrand as a function of A9o- The radial 
profile of the density is 



N f n f R f R 

Pn,i(r) = — J dA6 a J dr r f nA (r ) J dr r f nt i(r ) 

xSi(r' ,r ,A6 ) d Q e *(<t>-<t>'~A0 o ) ( 44 ) 
Jo 

The angular average of the last integral in the above for- 
mula has been worked out explicitly as elliptic integrals 
in the appendix, Eq. (B4). A few examples for vortex 
states in the lowest energy level, n = 0, are shown in 
Fig. ||. The density somewhat drops near the core, but 
generally speaking is quite uniform. The density never 
drops to zero at the center of the vortex. All the den- 
sity profiles depend logarithmically on the system size, 
and become completely uniform in the infinite system 
limit. This result is quite consistent with the conjecture 
that the density is roughly constant even inside the vor- 
tex core, proposed by FetterEil based on the argument of 
strong pair correlations in the wave function. In com- 
parisort-witb-,the Monte Carlo calculations for a static 
vortexJ13'tall3 my density profiles look like "spatially av- 
eraged", and miss several detailed features on the scale 
of particle spacing. For examples, the fractional density 
at the core center is over-estimated because my cut-off 
function g(r) is set to unity. Spatial oscillations due to 
higher-order correlation in the ground state are neglected 
by the Hartree approximation. The core size is also over- 
estimated possibly because the strong repulsion among 
particles has been under estimated in the dilute limit. In 
liquid helium, the range, in which the two-body poten- 
tial is strongly repulsive, is comparable to the spacing 
between two helium atoms. 

The current density in the azimuthal direction is 
jn,;(r), 




d 2 T, 



N 



ia_ 

rde 



(45) 



Once again, I can make an angular average so that the 
integrand depends only on A6q, 



n h 

A mr 



R 



dA6 / dr r f nt i(r ) / dr' Q r' f nA (r' ) 



R 



JO 

S, (r' , r , A0 Q )e- iAe ° J J dO e"^' ^ 



(46) 



The analytic form of the last angular integral has been 
work ed ou t in terms of elliptic integrals in the appendix, 
Eq. ( B16| ). A few plots of the azimuthal current density 
for vortex states in the lowest energy level are shown in 
Fig. |]. The current density is always zero at the origin 
as expected, and behaves like 1/r asymptotically. All 
current distributions are approximately independent of 
the system size. The vorticity is distributed around the 
vortex core over a region for about two or three atomic 
layers. 
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VI. LANDAU LEVELS AND EFFECTIVE MASS 

The motion of a vortex can now be understood by com- 
paring the weighting functions in Fig. |l| to the wave func- 
tions of an electron in an uniform magnetic field B. In 
the central gauge, A = (—yB/2,xB/2,0), the Hamilto- 
nian for the electron is 



H = 



(p-eA) 2 
2M 



P 

2M 



M /w\2 



(IP 



W T 

S 2 L *' 



(47) 



where u = \eB\/M is the cyclotron frequency, L z = 
xp y — yp x is the z component of the angular momen- 
tum operator, and s is the sign of eB. Without the 
L z term, the Hamiltonian represents a two-dimensional 
simple-harmonic oscillator, which has energy eigenvalues 
(n + l)huj/2, where n is a non-negative integer. The char- 
acteristic frequency for this oscillator is only half of the 
cyclotron frequency. The degeneracy for the n-th level is 
n+1. With the L z term, all energies with odd(even) n are 
shifted by an odd(even) multiples oifouj/2, and the energy 
spectrum changes to E = (n + l/2)hu. The degeneracies 
for the low-lying levels are now proportional to the sys- 
tem size. These energy levels are called Landau levels. In 
the polar coordinates, the radial wave function R nt m{r) 
in the n-th Landau level with angular momentum mh is 
a confluent hypergeometric function shown in Appendix 
A. There is also only one length scale set by t he streng th 
of the magnetic field. If this length scale, ^2h/\eB\, is 
chosen to be the unit length for the electron problem, the 
weighting function f n .N-m(fo) can be matched perfectly 
with the radial wave function Rn,sm(r). This one-to-one 
correspondence immediately tells us that the motion of 
a vortex is indeed cyclotron-like, and the degeneracy of 
each energy level in Fig. [2] is one per ir unit area, since 
the degeneracy of each Landau level is \eB\/h per unit 
area. 

Once the basis for the cyclotron motion of a vortex 
is established, the parameter, namely the effective mass, 
in the phenomenological theory can be derived from this 
microscopic theory. In cyclotron motion, the effective 
mass is related to the inverse of the energy-level spacing, 



M eff = 



2h 2 



Ae ~ a 2 Ae 



(48) 



where Ae is the energy-level spacing, and the role of eB 
is replaced by pn. p is the asymptotic mass density of the 
fluid, m/ira 2 , and k is the quantum of circulation, h/m. 
However, the energy spectrum in Fig. || is slightly devi- 
ated from the simple-harmonic spectrum that the level 
spacings tend to increase with energy. Each energy level, 
therefore, has a different effective mass. At present, it is 
not clear whether this deviation is due to the deficiency 
of the projection method, or is a real many-body effect. 
The single-particle picture is in fact broken by this small 
deviation because one cannot simply add an attractive 
potential in Eq. ( f47| ) to mimic the spectrum. The de- 
generacy of the Landau level will be broken in the same 



time by the additional potential. On the other hand, it 
is known, in the case of a nucleus, that the projection 
method in its simplest version only gives an approximate 
answer for. the mass. It was pointed out by Peierls and 
ThoulessB that the mass of a nucleus is only given ex- 
actly by working with further improved wave functions, 
which allow quantum fluctuations not only in the posi- 
tion but also in the velocity. In the case of a vortex, 
their suggestion is equivalent to taking the backfiow into 
account, and considering trial wave functions of the fol- 
lowing form, 



* = y d 2 rod 2 v G(ro,v )* r 



(49) 



The fluctuation in velocities is usually separable from 
the fluctuation in positions, and takes a simple Gaus- 
sian form. The effect of the backfiow is not clear at the 
moment, and is currently under investigation. 

The effective mass of a quantized vortex is a controver- 
sial subject even qualitatively.! 2 ^! 2 -!] One important issue 
is whether the mass scales with the system size. I have 
calculated the effective masses of the first few energy lev- 
els for various system sizes, and plot them in Fig. g. The 
data clearly shows that the effective mass scales loga- 
rithmically with the size of the vortex. This qualitative 
feature is not affected by the issue discussed in the end of 
the previous paragraph because the correction from the 
backfiow is clearly finite. This special scaling behavior is 
attributed to the long-ranged phase coherence in a super- 
fluid, which causes the logarithmic size-dependence of the 
overlap in Eqs. (36) to ( J3q ) between any two Feynman 
wave functions. It is not surprising that the overlap is 
significant non-zero only when the separation of two vor- 
tex center coordinates are close to each other. However, 
the functional form of the kernel for a short-ranged corre- 
lated system such as a nucleus is quite different because 
the width of the peak does not scale logarithmically with 
the system size. Even though the calculation is done rig- 
orously only for a system of weakly interacting bosons, 
I now argue that the scaling behavior is going to persist 
even in a strongly interacting system. The orthogonality 
is due to the phase factor associated with the vortex. Any 
two shifted vortex wave functions eventually lose their 
phase coherence at the large distance, even though they 
are only displaced by a small amount. Interaction may 
change the detail dependence on the vortex coordinates 
or the constants, but can only change the short-distance 
physics. For a sufficiently large vortex, the logarithm is 
going to be the dominate term, so is the scaling behav- 
ior. In a real system, the logarithm must be cut off at 
a certain length scale, which could be the distance that 
the condensate loses its phase coherence, or the spacing 
among vortices. 

In conclusion, I have presented a microscopic quantum 
theory to support the picture that a quantized vortex be- 
haves like an electron in a magnetic field. Vortex states 
with different angular momenta form highly degenerate 



8 



levels. The density is finite at the vortex core axis, and 
the vorticity is distributed inside the core. The dynami- 
cal effective mass defined as the inverse of the energy-level 
spacing scales logarithmically with the size of the vortex. 
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APPENDIX A: AN ELECTRON IN A MAGNETIC FIELD 



In this appendix I show the energy eigenfunctions of an electron in a uniform magnetic field.E3 I look for energy 
eigenfunctions in the polar coordinates, Rn,m(f)^ lme ■ The radial Schrodinger equation takes the following form, 



2M 



dr 2 



1 d 2ME (m Mujr\ 



dr h 2 \ r 2k 



■ 



Rn,m{r) = , 



(Al) 



where E = (n + l/2)fkj. It is obvious from the differential equation that the wave function falls off exponentially in 
the asymptotic region. In dimensionless variable £ = rj y / 2fi/\eB\, the radial wave function takes the following form, 

Rn,m(0 =9n,m(0e~ e/2 , (A2) 

where g n ,m(Ji) satisfies the differential equation, 



d 

- 2£ ) — + 2{2n + sm) - 



9n,m(0 = 



(A3) 



In looking for the series solutions of <7n,m(£), the indicial equation has two roots, ±to. Since the wave function is 
required to be regular at the origin, only the positive root satisfies the requirement, 



J=0 



where c, satisfies the following recursion relation, 



Cj+2 



\m\ +j — 2n— sm 
(\m\ +j + 2) 2 - m 2 



(A4) 



(A5) 



The energy is, therefore, quantized so that the infinite series terminates at the finite order j max = 2n + sm — |m|. 
This series can be rewritten in terms of the confluent hypergeometric function, 



9n,m(0 = £ |m| 1F1 [-(2n + sm - \m\)/2, \m\ + 1; £ 2 



where \Fi(a,b;x) is defined as 



iFi{a, b;x) = ^ 



(o + j -l)l x j 



(A6) 



(A7) 



Apart from the normalization constant, R n S m(0 is identical to the weighting function /„ at— m(^*o) i n Eq. (Q) as shown 
in Fig. 0. 



APPENDIX B: ANGULAR INTEGRALS 

This appendix contains several angular integrals for evaluating the constants in Eqs. (pi) and (|35|), and for calcu- 
lating the radial profile of the number density in Eq. ( ffi]) and the current density in Eq. (|f6|). The remaining radial 
parts of these integrals can then be carried out numerically in a reasonably length of time. 

The overlap integral in Eq. (^8|) can be written in the following form using complex variables, 



tfre* = / drr / dflj ^ = ^ffi^ = ^ , (Bl) 




^/^-vra^y- <B2) 



where 
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B = ^-e< C=-e l9 ° D=-± 
r' r 



„i9' n 



(B3) 



On the complex plane, there are two Riemann sheets and four branch points. The integration contour is along the 
unit circle, and there are always two branch points inside the contour and two points outside the contour. One can 
always rotate the axes, so that two branch points are lying on the s-axis. I choose the branch cuts in such a way 
that the contour is separated into two disjointed curves. The contour and the branch cuts are shown in Fig. ^. Using 
the formula for indefinite integrals obtained in the next appendix, I can carry out the angular part of the integral by 
suitably deforming the contours. Two additional contours are added to go around the branch points, and to avoid 
the branch cuts. The integral along the circular contour is now equal to the sum of four integrals on the straight 
segments plus the residue from a simple pole at the origin. The result is summarized as follows 



27r < /^e" 4Ae "/ 2 +4d 



I to 
>\ 

[r_o 
r, 



n(m,fe)-n(n 2 ,A;) 



2tt ^e- lAe °/ 2 - 4cJ II(ni, k) - U(n 2 , k) 



Vn 3 ' k 



for r > r ,r' Q 



1 1 



2tt, /4 e - jA <V 2 - 4C a U{n 3 ,k) - n(n 4 , k) 



Vn4 ' k 
for ro < r < r' 

for r < ro , r' 



(B4) 



where n(n, k) = T\(n,ir /2,k) is the complete elliptic integral defined in the next appendix. Various abbreviations in 
the equation are defined as follows 



ni 



n 2 



r 2 - r r' e~ iAe ° 

\/ {r 2 - rl){r 2 - r' 2 ) 

r' (r e~* AS °~r> ) 

r z — 7q 
r 2 (r -r' Q e tAe «) 

r o(r 2 - r' 2 ) 
-r 2 (r 2 + r' 2 - 2r a r' Q cos Afl ) 
(r 2 — r 2 )(r 2 — r' 2 ) 



(B5) 

(B6) 
(B7) 
(B8) 



and (113,124) are simply (714,712) with (ro,r' ) interchanging their roles. This is also the formula for calculating the 
number density in Eq. (|44|). 

The overlap integral with the kinetic energy term in Eq. (Bq) can be done in the similar way, 



drr / d9 



(re 16 - r Q e Wo )(re- 



»o) + (re 10 - r' Q e ie o)(re- lf > - r e- te °) 



2(re l6 — r Q e ie «)(re~ ie — r o e^ l0o )\re l6 — r e 1 ' 



018 



= I drr 



dz 
2^ 



IBC 



C 3 (z ~ A)(z — B) 
B(z — C) 3 (z — D) 3 



y/(z-A)(z-B)(z-C)(z-D) 
Using the same contour without the simple pole at the origin, I obtain the angular integral as 
4 



(B9) 
(BIO) 



< 



C 2 E(k) ~ iC 3 K(k) 
(1 - fc) 



^4G 



E(k)- 



Vk 



K(l/k) + VkE(l/k) 



for r>r ,r' Q r < r ,r' 
yk 

for ro < r < r' Q 



(Bll) 



where E(k) = E(tt/2, k) and K(k) = F(ir/2, k) are again the complete elliptic integrals defined in the next appendix, 
and 
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c 2 = 



y \r 2 - r 2 )(r 2 - r' 2 ) 
(r 2 - r r' e iAe o) 2 ' 
r r' sin A^o 



(r 2 - r r' Q e iAe ^)\J (r 2 - r 2 )(r 

To calculate the azimuthal current density in Eq. (|46|), I have the following integral, 

d 



-le 



-iA8 



f de 


P -i<f>' ( d jA 


' 2^ 


Km J 


dz 


1 z z 


Airz 


[z-A ' z-D 



00 



D 



C 



! (z-A)(z-B) 
(z-C){z-D) 



The result for r > ro , Tq is 



where 



ZI . /!^ e - 4 Ae /2 



C 4 



2C 4 



U(m,k) -U(n 2 ,k) +2C 5 E(k) 



r 2 - r r' Q e~ lAe ° 



r^J(r 2 -r 2 ){r 2 - r' 2 ) 

y \r 2 - r 2 )(r 2 - r' Q 2 ) 
r(r 2 — r r' Q e iAe<) ) 



(B12) 
(B13) 



(B14) 
(B15) 

(B16) 

(B17) 
(B18) 



For tq < r < r' and r < ro,r' , similar substitutions can be performed as in the previous cases, Eqs. (B4) and (Bll) 



APPENDIX C: FORMULA WITH ELLIPTIC INTEGRALS 

The following indefinite integrals are the basic ingredients in the previous appendix. These integral are calculated 
with the help of Mathematca. For the overlap integral in Eq. (B2), I use 



dz / (z — a)(z — b) 



2(a-b) 



z y (z - c)(z - d) y/( a -c)(b-d) 
where II(n, (f>, k) is the elliptic integral of the third kind, 



II(n, (j), k) — II(6n/a, 0, k) 



II(n, <p, k) — 



de 



o (1 - ri sin 2 6»)\/l - k 



sin 



and the parameters are defined as 



(CI) 

(C2) 

(C3) 
(C4) 

- (a-c)(b-d)- (C5) 
For the overlap integral weighted by the kinetic energy in Eq. ( B10| ), two indefinite integrals are needed. The first 



a — d 
b-d 



t (b-d)(z-a) 
\l(a-d)(z-b) 

(b-c)(a-d) 



one is 



dz 



\J \z — a)(z — b)(z — c)(z — d) \J (a — c)(b — d) 



F(^k) , 



(C6) 
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where F(<p, k) is the elliptic integral of the first kind 

Ffak) = 

The second one is 



(10 



o y/l — k sin 2 ( 



/ 



dz 



>{z-a){z-b) 



(z — c)(z — d) V (z — c){z — d) (c — d) 2 



(a-b)(c-d) 
y/(a-c)(b-d) 



F(cf>,k) 



-2^(a-c)(b-d)E( < f >) k) + (z- n} (z -c)(b-d) + (z~d)(b-c) 

\/{z — a)(z — b){z — c)(z — d) 

where E(<j>, k) is the elliptic integral of the second kind, 

E{<f>,k) = r d9\/l - k sin 2 e . 
Jo 



Finally, the integral for the distribution of the current density in Eq. ( B15 ) is 

z b c \ l(z-a)(z-b) _ 4y/(a-c){b-cC) 



dz 



\z — a 
4(a - b) 



z — d z — b z — c) \J (z — c)(z — d) 



c-d 



£7(0, A;) 



y/(a-c)(b-d) 



Il(n, <fi, k) — Tl(bn/a, 4>, k) 



2(z - a) [(z -c)(b-d) + {z- d)(b - c)] 
(c — d) \J{z — a)(z — b)(z — c)(z — d) 



(C7) 



(C8) 



(C9) 
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r [a] r [a] 

FIG. 1. The normalized weighting functions for angular momenta / = N, N — 1, N — 9, N — 10. These are the eigenfunctions 
of matrices with 250 x 250 elements at the system size R = 10. The normalization condition is J drf„ y i(r) 2 = 1. The solid 
lines are the corresponding radial wave functions of an electron. The only adjustable parameter for the solid curves is the 
normalization constant. 
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Log w (R/a) 

FIG. 2. The energy spectrum corresponding to the cy- 
clotron motion of a vortex relative to the lowest energy level. 
The dotted lines show the scaling of each energy level with 
the system radius R ranging from 20a to 10 a. The lowest 
energy level approximately scales as \xi{R/o) corresponding 
to the formation energy of a vortex. The same spectrum is 
obtained for any given angular momentum. 




r[o] 

FIG. 3. The radial profiles of the number density for three 
different states, each of which has total angular momentum 
N, N — 1, and N — 4, in the first Landau level. The curves 
represent the numerical results in a finite system with radius 
R — 10. The limit for all states in the infinite system is uni- 
form density. Instead of fixing the total number of particles, 
the normalization is chosen in such a way that the asymptotic 
value is unity. 



r[a] 

FIG. 4. The radial profiles of the azimuthal current den- 
sity for three different states, each of which has total angular 
momentum N, N — 1, and N — 4, in the first Landau level. 
The dashed lines represent the numerical results in a finite 
system with radius R — 10, and the solid lines represent the 
results for the infinite-system limit. The dotted line is the 
classical result for a singular vortex line. 




Log w (R/a) 

FIG. 5. The scalings of the inverse of energy-level spacings 
labeled as Ae„ — e n — £ n -i, which are related to the effective 
mass of a vortex. The symbols show the first five of them from 
top to bottom. The lines are the least-square fits, and clearly 
show that the effective mass scales logarithmically with the 
vortex size R. 
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FIG. 6. Branch points, A to D, branch cuts, AC and BD, 
and the integration contour deformed from the unit circle. 
This particular sequence of the branch points corresponds to 
the case r > ro , r' . There may be a simple pole located at 
the origin depending on which integrand is discussed. 
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